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Abstract 

A distributed single-hop wireless network with K links is considered, where the links are partitioned into a 
fixed number (M) of clusters each operating in a subchannel with bandwidth S. The subchannels are assumed 
to be orthogonal to each other. A general shadow-fading model, described by parameters [pt,w), is considered 
where a denotes the probability of shadowing and w (vj < 1) represents the average cross-link gains. The main 
goal of this paper is to find the maximum network throughput in the asymptotic regime of K — > oo, which is 
achieved by: i) proposing a distributed and non-iterative power allocation strategy, where the objective of each user 
is to maximize its best estimate (based on its local information, i.e., direct channel gain) of the average network 
throughput, and ii) choosing the optimum value for M . In the first part of the paper, the network throughput is 
defined as the average sum-rate of the network, which is shown to scale as 6(logA"). Moreover, it is proved 
that in the strong interference scenario, the optimum power allocation strategy for each user is a threshold-based 
on-off scheme. In the second part, the network throughput is defined as the guaranteed sum-rate, when the outage 
probability approaches zero. In this scenario, it is demonstrated that the on-off power allocation scheme maximizes 
the throughput, which scales as — log K. Moreover, the optimum spectrum sharing for maximizing the average 
sum-rate and the guaranteed sum-rate is achieved at M = 1. 
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I. Introduction 

A. History 

A primary challenge in wireless networks is to use available resources efficiently so that the network 
throughput is maximized. Throughput maximization in multi-user wireless networks has been addressed 
from different perspectives; resource allocation [3]-[5], scheduling [6], routing by using relay nodes [7], 
exploiting mobility of the nodes [8] and exploiting channel characteristics (e.g., power decay-versus- 
distance law [9]— [1 1], geometric pathloss and fading [12]— [14]). 

Among different resource allocation strategies, power and spectrum allocation have long been regarded 
as efficient tools to mitigate the interference and improve the network throughput. In recent years, 
power and spectrum allocation schemes have been extensively studied in cellular and multihop wireless 
networks [3], [4], [15]-[20]. In [19], the authors provide a comprehensive survey in the area of resource 
allocation, in particular in the context of spectrum assignment. Much of these works rely on centralized and 
cooperative algorithms. Clearly, centralized resource allocation schemes provide a significant improvement 
in the network throughput over decentralized (distributed) approaches. However, they require extensive 
knowledge of the network configuration. In particular, when the number of nodes is large, deploying such 
centralized schemes may not be practically feasible. Due to significant challenges in using centralized 
approaches, the attention of the researchers has been drawn to the decentralized resource allocation schemes 
[21]-[26]. 

In decentralized schemes, the decisions concerning network parameters (e.g., rate and/or power) are 
made by the individual nodes based on their local information. The local decision parameters that can be 
used for adjusting the rate are the Signal-to-Interference-plus-Noise Ratio (SINR) and the direct channel 
gain. Most of the works on decentralized throughput maximization target the SINR parameter by using 
iterative algorithms [23]-[25]. This leads to the use of game theory concepts [27] where the main challenge 
is the convergence issue. For instance, Etkin et al. [25] develop power and spectrum allocation strategies 
by using game theory. Under the assumptions of the omniscient nodes and strong interference, the authors 
show that Frequency-Division Multiplexing (FDM) is the optimal scheme in the sense of throughput 
maximization. They use an iterative algorithm that converges to the optimum power values. In [24], 
Huang et al. propose an iterative power control algorithm in an ad hoc wireless network, in which 
receivers broadcast adjacent channel gains and interference prices to optimize the network throughput. 
However, this algorithm incurs a great amount of overhead in large wireless networks. 

A more practical approach is to rely on the channel gains as local decision parameters and avoid 
iterative schemes. Motivated by this consideration, we study the throughput maximization of a distributed 
wireless network with K links, operating in a bandwidth of W. To mitigate the interference, the links 
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are partitioned into a fixed number (M) of clusters, each operating in a subchannel with bandwidth jfe, 
where the subchannels are orthogonal to each other. Throughput maximization of the underlying network 
is achieved by proposing a distributed and non-iterative power allocation strategy based on the direct 
channel gains, and then choosing the optimum value for M. 

B. Contributions and Relations to Previous Works 

In this paper, we study the throughput maximization of a spatially distributed wireless network with 
K links, where the sources and their corresponding destinations communicate directly with each other 
without using relay nodes. Wireless networks using unlicensed spectrum (e.g. Wi-Fi systems based on 
IEEE 802.11b standard [28]) are a typical example of such networks. The cross-link channel gains are 
assumed to be Rayleigh-distributed with shadow-fading, described by parameters (a,zu), where a denotes 
the probability of shadowing and w (w < 1) represents the statistical average of the Rayleigh distribution. 

The above configuration differs from the geometric models proposed in [8]-[ll], [29], in which the 
signal power decays based on the distance between nodes. Unlike [22]-[25] which relies on an iterative 
algorithm using SINR, we assume that each transmitter adjusts its power solely based on its direct channel 
gain. 

If each user maximizes its rate selfishly, the optimum power allocation strategy for all users is to 
transmit with full power. This strategy results in excessive interference, degrading the average network 
throughput. To prevent this undesirable effect, one should consider the negative impact of each user's 
power on other links. A reasonable approach for each user is to choose a non-iterative power allocation 
strategy to maximize its best local estimate of the network throughput. 

The network throughput in this paper is defined in two ways: i) average sum-rate and ii) guaranteed 
sum-rate. It is established that the average sum-rate in the network scales at most as Q(\ogK) in the 
asymptotic case of K — > oo. This order is achievable by the distributed threshold-based on-off scheme 
(i.e., links with a direct channel gain above certain threshold transmit at full power and the rest remain 
silent). Moreover, in the strong interference scenario, the on-off power allocation scheme is the optimal 
strategy. In addition, the on-off power allocation scheme is always optimal for maximizing the guaranteed 
sum-rate in the network, which is shown to scale as — log K. These results are different from the result 
of [30] where the authors use a similar on-off scheme for M = 1 and prove its optimality only among 
all on-off schemes. This work also differs from [31] and [32] in terms of the network model. We use 
a distributed power allocation strategy in a single-hop network, while [31] and [32] consider an ad hoc 
network model with random connections and relay nodes. 

We optimize the average network throughput in terms of the number of the clusters, M. It is proved 
that the maximum average sum-rate and the guaranteed sum-rate of the network for every value of a and 
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w is achieved at M = 1. In other words, splitting the bandwidth W into M orthogonal sub-channels does 
not increase the throughput. 

The rest of the paper is organized as follows. In Section HH the network model and objectives are 
described. The distributed on-off power allocation strategy and the network average sum-rate are presented 
in Section [TTTJ We analyze the network guaranteed sum-rate in Section [IV] Finally, in Section |Vj an 
overview of the results and some conclusion remarks are presented. 

C. Notations 

For any functions f{n) and g(n) [33]: 
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• f(n) pa g(n) means that f{n) is approximately equal to g{n), i.e., if we replace f{n) by g{n) in the 
equations, the results still hold. 

Throughout the paper, we use log(.) as the natural logarithm function and P{.} denotes the probability 
of the given event. Boldface letters denote vectors; and for a random variable x, x means E[x], where 
E[.] represents the expectation operator. RH(.) represents the right hand side of the equations. 

II. Network Model and Objectives 

A. Network Model 

In this work, we consider a single-hop wireless network consisting of K pairs of nodes[] indexed 
by {1,... operating in bandwidth W. All the nodes in the network are assumed to have a single 
antenna. The links are assumed to be randomly divided into M clusters denoted by C,, j = 1, ...,M 
such that the number of links in all clusters are the same. Without loss of generality, we assume that 
Cj = {(j — l)n + 1, jn}, where n = j| denotes the cardinality of the set Cj which is assumed to be 
known to all userj^. To eliminate the mutual interference among the clusters, we assume an M -dimensional 
orthogonal coordinate system in which the bandwidth W is split into M disjoint subchannels each with 

'The term "pair" is used to describe a transmitter and its corresponding receiver, while the term "user" is used only for the transmitter. 
2 It is assumed that K is divisible by M, and hence, n = p- is an integer number. 
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bandwidth § . It is assumed that the links in Cj operate in subchannel j. We also assume that M is fixed, 
i.e., it does not scale with K. The power of Additive White Gaussian Noise (AWGN) at each receiver is 
^r¥-, where N is the noise power spectral density. 

The channel model is assumed to be flat Rayleigh fading with the shadowing effect. The channel gaii|] 
between transmitter k and receiver i is represented by the random variable Cki- For k = i, the direct 
channel gain is defined as Cu — ha where ha is exponentially distributed with unit mean (and unit 
variance). For k ^ i, the cross channel gains are defined based on a shadowing model as follows: 

r A J Pkihki, with probability a 

I 0, with probability 1 — a, 
where h^s have the same distribution as ha's, < a < 1 is a fixed parameter, and the random variable 
Pki, referred to as the shadowing factor, is independent of and satisfies the following conditions: 

• A™ < Pki < Anax, where (3 min > and /3 max is finite, 
. E[&J ±w< 1. 

It is also assumed that {£fcj} and {Pki} are mutually independent random variables for different (k,i). 

All the channels in the network are assumed to be quasi-static block fading, i.e., the channel gains 
remain constant during one block and change independently from block to block. In addition, we assume 
that each transmitter knows its direct channel gain. 

We assume a homogeneous network in the sense that all the links have the same configuration and use 
the same protocol. We denote the transmit power of user i by pi, where pi E 2? = [0,P maa ]. The vector 
pC?) — {jp^__^ n+l ^ ...,pj n ) represents the power vector of the users in Cj. Also, P_] denotes the vector 
consisting of elements of P^ other than the i th element, i G Cj. To simplify the notations, we assume 
that the noise power ^jf- is normalized by V max . Therefore, without loss of generality, we assume that 
P max = 1. Assuming that the transmitted signals are Gaussian, the interference term seen by link i G Cj 
will be Gaussian with power 

h = ^ CkiPk- (2) 

fceCj 

k^i 

Due to the orthogonality of the allocated sub-channels, no interference is imposed from links in Ck on 
links in Cj, k ^ j. Under these assumptions, the achievable data rate of each link i G Cj is expressed as 

/^UF>) = |log(l + I ^), (3) 

where Cf* = (C((j-i) n+ i)i, ...,£(j n )i). To analyze the performance of the underlying network, we use the 
following performance metrics: 

3 In this paper, channel gain is defined as the square magnitude of the channel coefficient. 
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Network Average Sum-Rate: 



p A iu< 



M 



j=i leCj 



(4) 



where the expectation is computed with respect to C\ j \ This metric is used when there is no decoding 
delay constraint, i.e., decoding is performed over arbitrarily large number of blocks. 
• Network Guaranteed Sum-Rate: 



M 



j=i leCj 



(5) 



in which for all h tt , I G C 7 , we have 



R*(h u ) = su V R(hu) 



such that 



p{i?KP 0) ,A a) )<^)}^o. 



(6) 



(7) 



This metric is useful when there exists a stringent decoding delay constraint, i.e, decoding must be 
performed over each separate block, and a single-layer code is used. In this case, as the transmitter 
does not have any information about the interference term, an outage event may occur. Network 
guaranteed throughput is the average sum-rate of the network which is guaranteed for all channel 
realizations. 



B. Objectives 

Part I: Maximizing the network average sum-rate: The main objective of the first part of this paper 
is to maximize the network average sum-rate when the interference is strong enough, i.e., E[ij] = 
This is achieved by: 

- Proposing a distributed and non-iterative power allocation strategy, where each user maximizes its best 
estimate (based on its local information, i.e., direct channel gain) of the average network sum-rate. 

- Choosing the optimum value for M. 

To address this problem, we first define a utility function for link i e Cj (j — 1, M) that describes 
the average sum-rate of the links in cluster Cj as follows 



]>>(po-UP) 



leCi 



(8) 



where the expectation is computed with respect to {Cki}k,ieCj excluding k — I — i (namely ha). As 
mentioned earlier, ha is considered as the local (known) information for link i, however, all the other 
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gains are unknown to user i which is the reason behind statistical averaging over these parameters in ©. 
User % selects its power using 



pi = arg max Ui(p h ha 



(9) 



It will be shown that when the number of links is large and the interference is strong enough, the optimum 
power allocation strategy for the optimization problem in © is the on-off power scheme. Assuming that 
the channel gains change independently from block to block, each user updates its on-off decision based 
on its direct channel gain in each block. Given the optimum power vector P = (p(j-i) n +i, ■■■,Pjn) 
obtained from ©, the network average sum-rate is then computed as ©. Next, we choose the optimum 
value of M such that the network average sum-rate is maximized, i.e., 



M 



arg max R a 

M 



(10) 



Also, for the moderate and the weak interference regimes (i.e., E[ij] = 0(1)), we obtain upper bounds 
for the network average sum-rate. 

Part II: Maximizing the network guaranteed sum-rate: The main objective of the second part is 
finding the maximum achievable network guaranteed sum-rate in the asymptotic case of K — > oo. For 
this purpose, a lower bound and an upper-bound on the network guaranteed sum-rate are presented and 
shown to converge to each other as K — > oo. Also, the optimum value of M is obtained. 



III. Network Average Sum-rate 
A. Strong Interference Scenario (E[ij] = co(l)) 

In order to maximize the average sum-rate of the network, we first find the optimum power allocation 



policy. Using ([8]), we can express the utility function of link % £ Cj, j = 1, M, as 

i(pi, ha) = Ri(pi, ha) + ^ Ri(Pi), 



U; 
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Ri(Pi, ha) = E 
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with the expectation computed with respect to ij defined in ©, and 
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log 1 + 
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(12) 

(13) 
(14) 

(15) 



with the expectation is computed with respect to pff] and {C k i} kleC excluding I = It is worth 
mentioning that the power p ri in (TT3T ) prevents the ith user from selfishly maximizing its average rate 
given in (fT2l) . Using the fact that all users follow the same power allocation policy, and since the channel 
gains Cki are random variables with the same distributions, Ri(pi) becomes independent of I. Thus, by 
dropping the index I from R~i(pi), the utility function of link i can be simplified as 

Ui(pi, ha) = Ri(pi, ha) + (n - l)R{pi). (16) 

Noting that pi depends only on the channel gain ha, in the sequel we use pi = g(hu). 

Lemma 1 Let us assume ~&[pk] — q n > < a < 1 is fixed and the interference is strong enough (E[Jj] = 
u){\)). Then with probability one (w. p. 1), we have 

li ~ (n - l)aq n , (17) 

as K — » oo (or equivalently, n — > oo), where a = aw. More precisely, substituting li by (n — l)aq n does 
not change the asymptotic average sum-rate of the network. 

Proof: See Appendix HI ■ 

Lemma 2 For large values of n, the links with a direct channel gain above hxh — c log n, where c > 1 
is a constant, have negligible contribution in the network average sum-rate. 

Proof: See Appendix HO ■ 
From Lemma [2] and for large values of n, we can limit our attention to a subset of links for which the 
direct channel gain ha is less than clogn, c > 1. 

Theorem 1 Assuming the strong interference scenario and sufficiently large K, the optimum power 
allocation policy for ((9]) is pi = g{hu) = U{hu — r n ), where r n > is a threshold level which is a 
function of n, and U(.) is the unit step function. Also, the maximum network average sum-rate in (0) is 
achieved at M = 1 and is given by 

W 

Rave^— log K. (18) 

a 

Proof: The steps of the proof are as follows: First, we derive an upper bound on the utility function 
given in (fT6l) . Then, we prove that the optimum power allocation strategy that maximizes this upper 
bound is pi = g(hu) = U(hu — r„). Based on this power allocation policy, in Lemma 0] we derive the 
optimum threshold level r n . We then show that using this optimum threshold value, the maximum value 

4 Note that the power of the users are random variables, since they are a deterministic function of their corresponding direct channel gains, 
which are random variables. 
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of the utility function in (fT6b becomes asymptotically the same as the maximum value of the upper bound 
obtained in the first step. Finally, the proof of the theorem is completed by showing that the maximum 
network average sum-rate is achieved at M — 1. 
Step 1: Upper Bound on the Utility Function 



Let us assume E \p k ] = q n . Using the results of Lemma [H Ri(pi, ha) in (fT6l) can be expressed as 
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(21) 



In the above equations, (a) follows from the fact that ha is a known parameter for user % and pi = g(ha) 
is the optimization parameter. With a similar argument, (TT5T ) can be simplified as 
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(23) 
(24) 



as K — > oo, where the expectation is computed with respect to ft^, p/ and 0u, and A' = (n — 2)aq n + 
Also, (a) comes from the shadowing model described in (OQ). Using (1201) . (l24l) . and the inequality 



AW 

n 

log(l + x) < x, the utility function in (fT6l) is upper bounded asj 

hupi 
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Noting that h u is independent of h u , i I, we have 
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where 

(jl = E [hupi] , 

5 Note that the factor (n — 1) in dl6t is replaced by n in l !25t . which does not affect the validity of the equation. 



(25) 

(26) 
(27) 
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and Ei(x) = — f_ ^-dt, x < is the exponential-integral function [34]. Thus, the right hand side of 
(|25l ) is simplified as 

W [i 



Ui{Pi, ha) < Jtf-^Pi - n ^f E 



1 -^-^ ( A' 

-ePuPi Ei 



+ nil — a) 



MX' 



(29) 



fiilVi \ PilPi 

where the expectation is computed with respect to (3u . An asymptotic expansion of Ei(x) can be obtained 
as [34, p. 951] 
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as x — > — oc. Setting L = 4, we can rewrite (1291) as 
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the term E 
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d 2 Ei(pi,hii) 



n- 



dp. 



Step 2: Optimum Power Allocation Policy for Ej(pj, ha) 

Using the fact that G [0,1], the second-order derivative of (|32|) in terms of pi, 

< ^tt~ f t^- -rlrPi I » 1S positive^ as A' — > oo. Thus, ((321) is a convex function of pj. It is known that a 

MA' \X' 2 X' 6 J 

convex function attains its maximum at one of its extreme points^ of its domain [35]. In other words, the 
optimum power that maximizes (|32l) is pi E {0, 1}. To show that this optimum power is in the form of a 
unit step function, it is sufficient to prove that pi = g(hu) is a monotonically increasing function of ha. 

Suppose that the optimum power that maximizes Sj(pj, ha) is pi = 1. Also, let us define h u = ha + 6, 
where 5 > 0. From (1321) . it is clear that Ej(pi, hu) is a monotonically increasing function of ha, i.e., 

Ei(pi = 1, /4) > = 1, /iii). (34) 
On the other hand, since the optimum power is p- L = 1, we conclude that 

Ei(pi = 1, > Ei(pi = 0, /iii). (35) 
Using the fact that Hj(pj = 0, = S^Pi = 0, h'^), we arrive at the following inequality 

Ei(pi = 1, h' u ) > Ei(pi = 0, h' u ). (36) 

6 It is observed from OOb and Q2b that for any value of L > 4, the second-order derivative of 02b in terms of pt is positive too. 
7 In the power domain 3 s = [0, 1], the extreme points are and 1. 
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From (|34l)-(l36lh it is concluded that g(hu) is a monotonically increasing function of hu. Consequently, 
the optimum power allocation strategy that maximizes Eifa, hu) is a unit step function, i.e., 

A J 1, if ha > r n 

Pi= < (37) 
I 0, Otherwise, 

where r n is a threshold level to be determined. We call this the threshold-based on-off power allocation 
strategy. It is observed that the optimum power pi is a Bernoulli random variable with parameter q n , i.e., 

m) = I qn) pl = h 08) 

{ 1 - q n , pi = 0, 

where /(.) is the probability mass function (pmf) of pi. We conclude from (|37T) and (1381) that the probability 
of link activation in each cluster is q n = P {hu > r n } = e~ Tn which is a function of n. 
Step 3: Optimum Threshold Level r n 

From Step 1, it is observed that for every value of pi we have 

Ui{pi,hu) < Ei(pi,hu). (39) 

The above inequality is also valid for the optimum power pi obtained in Step 2. Thus, using the fact that 
for X < Y, E[X] < E[Y], we conclude 

E[ui{pi,hii)} < E[Ei{pi,hii)}, (40) 

where the expectations are computed with respect to ha. In the following lemmas, we first derive the 
optimum threshold level r n that maximizes E[Hj(pj, ha)}, and then prove that this quantity is asymptotically 
the same as the optimum threshold level maximizing^ E[ui(pi, hu)], assuming an on-off power scheme. 
We also show that the maximum value of E[ui(pi,ha)] (assuming an on-off power scheme) is the same 
as the optimum value of E[S;(p;, hu)], proving the desired result. 

Lemma 3 For large values of n and given < a < 1, the optimum threshold level that maximizes 
E[Hj(pj, ha)] is computed as 

f n ~logn. (41) 
W 

Also, the maximum value of E[Hj(pi, hu)] scales as -^rlogn. 

Proof: See Appendix [Till ■ 

Lemma 4 For large values of n and given < a < 1, 

8 In fact, since the threshold r n is fixed and does not depend on a specific realization of hu, finding the optimum value of r n requires 
averaging the utility function over all realizations of ha. 
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i) The optimum threshold level that maximizes E[ui(pi, ha)) is computed as 

f n = logn-21oglogn + 0(l), (42) 

ii) The probability of link activation in each cluster is given by 

q n = 5 , (43) 

n 

where 5 > is a constant, 

W 

Hi) The maximum value ofE[ui(pi, ha)] scales as -j^z logn. 

Proof: See Appendix HVl ■ 
Step 4: Optimum Power Allocation Strategy that Maximize Ui(pi, ha) 

In order to prove that the utility function in (fT6l) is asymptotically the same as the upper bound Sj(j?j, ha) 



obtained in (|32l) . it is sufficient to show that the low SINR conditions in (|20l) and (|24j) are satisfied. Using 

([20]). (EB and (03]), the SINR is equal to where 

A 

n N Q W 

A ^ aS \og z n + — ■ (44) 

It is observed that A goes to infinity as n — > oo. On the other hand, since we are limiting our attention 
to links with ha < h Th = clogn, we have 

huPi ofJ-V (45) 



A \ log n J 

when n —>■ oo. Thus, for large values of n, the low SINR condition, ^p- <C 1, is satisfied. With a 
similar argument, the low SINR condition for (|24j) is satisfied. Hence, we can use the approximation 
log(l + x) ~ x, for x <C 1, to simplify (|20l) and (1241) as follows: 



Ri{Pi,hii) « j^yPi, (46) 



aW 7 



+ (47) 



.PiihuPi + A' 

Consequently, the utility function is the same as the upper bound Sj(j)j, /ijj) obtained in (|32l) . 

when n — > oo. Thus, the optimum power allocation strategy for © is the same as the optimum power 
allocation policy that maximizes ^(pi, ha). 
Step 5: Maximum Average Network Sum-rate 

Using ®, the average utility function of each user i, E [uiip^ ha)] , i £ C,-, is the same as the average 
sum-rate of the links in cluster Cj represented by 

flffi = E E , j = 1, M. (48) 
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- U) 

where P is the on-off powers vector of the links in cluster Cj. In this case, the network average sum-rate 
defined in © can be written as 

M 

Rave = ^ ^ Rave (49) 

i=i 

(a) Wf n 



a 



(50) 



K 

where (a) follows from (1D-20I) of Appendix ITVl Using (|42l) . and noting that n = — , we have 



Step 6: Optimum Spectrum Allocation 



W 1 K 

Rave ~—\0g—. (51) 

a M 



According to (1501) . the network average sum-rate is a monotonically increasing function of f n . Rewriting 
equation (ID- 151) of Appendix [TV] which gives the optimum threshold value for the on-off scheme: 

- e*> log ( 1 + + l + = 0, (52) 

\ na J na + r n e Tn 

it can be shown thaj^l 

fle T ' n w not, (53) 

which implies that f n is an increasing function of n. Therefore, the average sum-rate of the network is an 
increasing function of n and consequently, noting that n = is a decreasing function of M. Hence, the 
maximum average sum-rate of the network for the strong interference scenario and < a < 1 is obtained 
at M = 1 and this completes the proof of the theorem. ■ 
Motivated by Theorem [Q we describe the proposed threshold-based on-off power allocation strategy 
for single-hop wireless networks. Based on this scheme, all users perform the following steps during each 
block: 

1- Based on the direct channel gain, the transmission policy is 

{I, if ha > r n 
0, Otherwise. 

2- Knowing its corresponding direct channel gain, each active user i transmits with full power and rate 

* = b s( 1+ („-i )a X + ^ )- (54) 

3- Decoding is performed over sufficiently large number of blocks, yielding the average rate of ■J^ log K 
for each user, and the average sum-rate of ¥ log K in the network. 

'in deriving l !53t . we have used the fact that T " e . " <C 1, which is feasible based on the solution given in {42}. 
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Remark 1- Theorem Q] states that the average sum-rate of the network for fixed M depends on the 
value of a = aw and scales as log j|. Also, for values of M such that log M = o(log K), the network 
average sum-rate scales as ¥ log K. 

Remark 2- Let rrij denote the number of active links in Cj. Lemma 4 states that the optimum selection 
of the threshold value yields E[rrij] = nq n = (log 2 n) . More precisely, it can be shown that the optimum 
number of active users scales as 6 (log 2 n) , with probability one. 



B. Moderate and Weak Interference Scenarios (E[lj] = 0(1)) 

Theorem 2 Let us assume K is large and M is fixed. Then, 

i) For the moderate interference (i.e., E[ij] = 0(1)), the network average sum-rate is bounded by 
Rave < O(logn). 

ii) For the weak interference (i.e., E[Ij] = o(l)), the network average sum-rate is bounded by R ave < 
o(logn). 



Proof: i) From ©, we have 

Rn.iw. 



w 

M 



log 1 



(a) 
< 



< 



(b) 
< 



log 1 + 



M 

j=i iec 3 

> > — E 

M 

M TI/ 

x - x - W ' 

> > — E 

j=i leCj 

EE -login 



3=1 l&Cj 



M 




log 1 + 



M cM 

< irr n( ln log n, 
iv 



(55) 

(56) 

(57) 

(58) 
(59) 



where (a) follows from Lemma 2, which implies that the realizations in which hu > clogn for some 
c > 1 has negligible contribution in the network average sum-rate, (b) results from the Jensen 's inequality, 
E [logx] < log(E [x]), x > 0. Also, (c) follows from the fact that log(l + x) < x, x > 0. Since for the 
moderate interference, E[lj] = anq n = 0(1), and using the fact that M is fixed, we come up with the 
following inequality 

cM 



Rave — 



aN Q 
6(logn 



9(1) logn 



(60) 
(61) 
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ii) For the weak interference scenario, where E[ij] = anq n = o(l), and similar to the part (i), it is 
concluded from (|59l ) that 

Rave < 4rro(l)logn (62) 

= o(logn). (63) 

■ 

Remark 3- It is concluded from Theorems Q] and [2] that the maximum average sum-rate of the proposed 
network is scaled as B(\ogK). 

C. M Not Fixed (Scaling With K) 

So far, we assume that M is fixed, i.e., it does not scale with K. In the following, we present some 
results for the case that M scales with AO It should be noted that the results for M = o(K) is the same 
as the results in Theorem \T\ 

Theorem 3 In the network with the on-off power allocation strategy, if M = <d(K) and < a < 1, then 
the maximum network average sum-rate in @ is less than that of M = 1. Consequently, the maximum 
average sum-rate of the network for every value of 1 < M < K is achieved at M = 1. 

Proof: See Appendix IVl ■ 
Remark 4- According to the shadow-fading model proposed in (OQ), it is seen that for a = 0, with 
probability one, Cki = 0, k i. This implies that no interference exists in each cluster. In this case, the 
maximum average sum-rate of the network is clearly achieved by all users in the network transmitting at 
full power. It can be shown that for every value of 1 < M < K, the maximum network average sum-rate 
for a = is achieved at M — 1 (See Appendix [VTl for the proof). 

Remark 5- Noting that for M = K only one user exists in each cluster, all the users can communicate 
using an interference free channel. It can be shown that for M = K and every value of < a < 1, the 
network average sum-rate is asymptotically obtained as 

R ave » W{\og K - log N W - 7 ), (64) 

where 7 is Euler's constant (See Appendix IVIII for the proof). Therefore, for every value of < a < 1, it 



is observed that the average sum-rate of the network in (|64l) is less than that of M = 1 obtained in (|T8l) . 

Remark 6- Note that for M = 1, in which the average number of active links scales as 6 (log 2 K) (in 
the optimum on-off scheme), we have significant energy saving in the network as compared to the case 
of M = K, in which all the users transmit with full power. 

'"Obviously, we consider the values of M which are in the interval [1, K\. 
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D. Numerical Results 

So far, we have analyzed the average sum-rate of the network in terms of M and a, in the asymptotic 
case of K — > oo. For finite number of users, we have evaluated the network average sum-rate versus the 
number of clusters (M) through simulation. For this case, we assume that all the users in the network 
follow the threshold-based on-off power allocation policy, using the optimum threshold value. In addition, 
the shadowing effect is assumed to be lognormal distributed with mean va < 1 and variance 1. Fig. \T\ 
shows the average sum-rate of the network versus M for K = 20 and K = 40, and different values of 
a and zu. It is observed from this figure that the average sum-rate of the network is a monotonically 
decreasing function of M for every value of (a,w), which implies that the maximum value of R ave is 
achieved at M = 1. 

Based on the above arguments, we have plotted the average sum-rate of the network versus K for 
M — 1 and different values of (a,w). It is observed from Fig. [2] that the network average sum-rate 
depends strongly on the values of (a,zu). 

IV. Network Guaranteed Sum-Rate 

Recalling the definition of the network guaranteed sum-rate in ©, in this section we aim to find the 
maximum achievable guaranteed sum-rate of the network, as well as the optimum power allocation scheme 
and the optimum value of M. 

Theorem 4 The guaranteed sum-rate of the underlying network in the asymptotic case of K — > oo is 
obtained by 

W 

R g ~—\ogK, (65) 

a 

which is achievable by the decentralized on-off power allocation scheme. 

Proof: In order to compute the guaranteed rate for link I £ Cj, we first define the corresponding 
outage event as follows: 



In the following, we give an upper-bound and a lower-bound for R g and show that these bounds converge 
to each other as K — >• oo (or equivalently, n — » oo). 




(66) 



(67) 
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(b) 

Fig. 1. Network average sum-rate vs. M for a) K — 20, a — 1, 0.5, 0.1 and shadowing model with za 
and b) K = 40, a — 0.5 and shadowing model with xn = 1, 0.4, 0.1 and variance 1. 



= 0.5 and variance 1, 
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30 




(b) 

Fig. 2. Network average sum-rate vs. K for M = 1 and a) shadowing model with za = 0.5 and variance 1, and a = 1, 0.7, 0.4, 0.1, 
and b) shadowing model with w = 1, 0.7, 0.4, 0.1 and variance 1, and a = 0.5. 
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Upper-bound: An upper-bound on the guaranteed sum-rate can be given by lower-bounding the outage 
probability as follows: 



P 



> p 



M 



in which we have used the fact that log(l + x) < x. Denoting v = tin, we can write 



(?) 



(a) 

> P<| r 

(b) 

> 1 



-mu)B.(u) < e e^)(^flM-Pi«')' 



(68) 
(69) 

(70) 
(71) 



for some positive In the above equation, (a) results from (l69l) , noting that £(z/) > 0, and (6) follows 
from Markov's inequality [36, p. 77], and the expectation is taken with respect to l\. The above equation 
implies that finding an upper-bound for E [e~ Il ^ R ^~\ is sufficient for the lower-bounding the outage 
probability. For this purpose, using ©, we can write 



E [e 



-hi{u)R{v) 



E 



(a) 



(6) 



(c) 



ri E [< 



keCj 

k^l 



ri E [< 



keCj 



n-1 



k ^ I. 



(72) 
(73) 

(74) 

(75) 



In the above equation, (a) follows from the fact that {£ki}kec with k ^ I, and {pfcj-fcec, are mutually 
independent random variables, (6) results from writing Cm as uufikihki (from (OQ)), in which is an 
indicator variable which takes zero when Cm = and one, otherwise, (c) follows from the symmetry 
which incurs that all the terms E \e-^ R{ -^ Uklf3klhklPk '\ , k 6 C,, are equal. Noting that u kU j3 k u h k u and 
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p k are independent of each other, we have 



E [e 



-i{v)R{v)u k if3 k ih k ip k 



] = E/»„ [®h kl [E^ [E Pk [e 



-£,{v)R(v)u k if3 k ih kl p k 



]]]] 



(a) 



< %„ [(1 - ?n) + q n e 



]]] 



(6) 



E Pkl [E hkl [(1 - g n ) + q n (l - a + ae 



-t{v)R{u)l3 kl h kl 



)]] 



= Eft, 
< 1 - 

(e) 



1 - ag n + 
1 - 



1 + w £(i/)fl(i/) 
otq n f3kii{v)R{v) 



1 + kl £(u)R(p) 
aq n wi(p)R(v) 
l + (3 m ^{u)R{u) 

■■qnt(v)R(.v) 



(76) 
(77) 
(78) 
(79) 

(80) 

(81) 

(82) 



In the above equation, (a) follows from the fact that e~ e ' x < (1 — x) + xe~ e , V6> > and < x < 1, 
noting that E[p fc ] = g n . (6) results from the definition of u k i, which is an indicator variable taking zero 
with probability 1 — a and one, with probability a. (c) follows from the fact that as hu is exponentially- 
distributed, we have E hkl [ e ~^ )RMl3klhkl ] = l+j3kl ^ v)R{y) ■ (d) results from the facts that (3 M < /3 max and 
E[Pki) — tt. Finally, (e) follows from the fact that 1 — x < e~ x , Vx and noting that aw = a. 
Combining (1751) and (f8~2l) and substituting into (1711) yields 



(83) 
(84) 



where t(u) = (n _ 1)6 ^ 

Consider the cases of E{//} = u;(l) (strong interference) or E{^} = 0(1) (moderate interference). 
Let us define 7 4 mm ( 1} M-^). Setting £(,) 4 ^ ^^J we have TT £f^ + 

M = (n - l)ag„ + (1 - |)^, and as a result, 



p{Op)} > l- e /W^-l)**,-^] ^ « 



> 1— e 2M /3maxl KM/. 

Since = 6(1), it follows that the necessary condition to have 

Z1V1 Pmax 

t(y) = - — — Pl " ^. NqW . In other words, 



R*(») < 



Pi" 



(85) 
(86) 

is having R(u) < 



(87) 
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which implies that R g defined in © is upper bounded by 



R g < nWE, 



(n-l)ag n + (l-|)^ 
nWE v \piu] 



(n-l)ag n + (l-i)^- 



Now, defining \P n = log n + 2 log log n, we have 



E„M < E\p l v\v<V n ]W{v<V n } + -E\p l u\v>* n ] 

(a) 

< q n y n + E[u\u>V n ]F{u>y n } 

i gA + (^+l)e"*" 

to 



(88) 
(89) 

(90) 
(91) 
(92) 
(93) 



g n logn. 

In the above equation, (a) comes from the facts that 

E \jnu\u < #„] F{u < * n } < V n E [ Pl \u < * ft ] P{// < * n } < * n E[p z ] = V n q n , 

and < p; < 1. (6) results from the fact that v is exponentially-distributed, (c) follows from the facts that 
i) as we are considering the strong and moderate interference scenarios, it yields that (n — l)aq n = ^(1), 
or equivalently, q n = fi(^), and ii) the term (\l/ n + l)e~^ n scales as ralc | gw (due to the definition of \l/ n ) 
which is negligible with respect to the first term q n ^ n . Combining (|89l) and (|93l yields 

Wnq n log n 



^ < 



(n - l)o;g n + (1 



7\ N W 
2 / M 



w 

< —\ogn 

a 

W 

a 



In the case of weak interference, we have 



E\piv] 



M 



Mn 



E[p,z/]. 



Rewriting (1921) . we obtain 



EM < ?A + (*»+l)e A , V^ n >0. 



Selecting ^ n = log(g n 2 ) and defining e = nq n , we have 



^ ^ '^T ( l °Sn~\og(e x )) . 

iVn 



(94) 

(95) 
(96) 

(97) 
(98) 

(99) 

(100) 
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As in the weak interference scenario we have e = o(l), it follows from the above equation that R g 
o(W\ogn) in this scenario. Comparing with (|96l ), it follows that 

W 



R 9 <—\ogK. 

a 



(101) 



Lower-bound For the lower-bound, we consider the on-off power allocation scheme with r n = log n — 
2 log log n. Also, assume that M = 1 (or equivalently, n = K). Noting q n = e~ Tn , we obtain 

E[J,] = {n- l)aq n = e(\og 2 n). (102) 

Therefore, using the result of Lemma 1, it is realized that with probability one (n — l)aq n (l — e) < Ii < 
(n — l)aq n (l + e), for some e = o(l). In other words, defining 



$(v) = log 1 + 



in - l)ag n (l + e ) + M 
it follows that 

which implies that R*(v) > <&(u). As a result, 
i? g > nWE[<P(i/)] 



o(l), 



nWE 



log 1 + 



(n-l)ag n (l + e) + ™ 



(a) 



^/ log 1 + 

\ (n-l)ag„(l + e 

> nW I log I 1 + 



M 



V 



n-l)aq n (l + e) + ^ 



e v dv 



e v du, 



(103) 
(104) 

(105) 
(106) 

(107) 

(108) 



M 



where \l/ n = log n + 2 log log rt and (a) follows from the on-off power allocation assumption. As [n 



l)aq n (l + e) = 6 (log 2 n), it follows that 



that 



log 1 + 



[n - l)aq n (l + e) + 



o(l) in the interval [r n , \l/ n ], which implies 



(109) 



(n-l)afc(l + e) + M 
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in the interval of integration [r n , \& n ]. Hence, 



— 7 ; ; ; m~w I ve dv (HI) 

(n-l)a ?n (l + e) + Mi rn 

= 7 ir7T —- W ((^ + l)e--(*n + l)e-*") (112) 



M 

(a) nWr n q n 

{n - l)aq n {l + e) + « 



(113) 



M 

~ — logn (114) 

a 

W 

= —logK, (115) 

a 

where (a) results from the facts that (\P n + l)e~*™ <C (r n + l)e~ Tn and e~ Tn = g n . Combining the above 
equation with (11011) . the proof of Theorem 0] follows. ■ 
Remark 7- Similar to the proof steps of Theorem 1, it can be shown that the optimum value of 
M is equal to one. In fact, since the maximum guaranteed sum-rate of the network is achieved in the 
strong interference scenario in which the interference term scales as naq n with probability one, it follows 
that the maximum network average sum-rate and the network guaranteed sum-rate are equal. Therefore, 
the optimum spectrum sharing for maximizing the network guaranteed sum-rate is the same as the one 
maximizing the average sum-rate of the network (M = 1). 

V. Conclusion 

In this paper, a distributed single-hop wireless network with K links was considered, where the links 
were partitioned into a fixed number (M) of clusters each operating in a subchannel with bandwidth 
The subchannels were assumed to be orthogonal to each other. A general shadow-fading model, described 
by parameters (a, vj), was considered where a denotes the probability of shadowing and w (vj < 1) 
represents the average cross-link gains. The maximum achievable network throughput was studied in the 
asymptotic regime of K — » oo. In the first part of the paper, the network throughput is defined as the 
average sum-rate of the network, which is shown to scale as 0(logA'). Moreover, it was proved that in 
the strong interference scenario, the optimum power allocation strategy for each user was a threshold- 
based on-off scheme. In the second part, the network throughput is defined as the guaranteed sum-rate, 
when the outage probability approaches zero. In this scenario, it was demonstrated that the on-off power 
allocation scheme maximizes the network guaranteed sum-rate, which scales as ^- log K. Moreover, the 
optimum spectrum sharing for maximizing the average sum-rate and guaranteed sum-rate is achieved at 
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Appendix I 
Proof of Lemma Q] 

Let us define Xk — £kiPk, where Cki is independent of pk, for k i. Under a quasi-static Rayleigh 
fading channel model, it is concluded that Xk s are independent and identically distributed (i.i.d.) random 
variables with 



Var [xk] 



E [CkiPk] = aq n , 
E [x\] - E 2 [xk] 



(a) 

< 2anq n — (aq. 



2 

n) ; 



(A-l) 
(A-2) 
(A-3) 



where E [h^] = 2 and a = aw. Also, (a) follows from the fact that p\ < p/.. Thus, E[p|] < E[pfc] = g n . 



The interference U = ^keCj Xk is a random variable with mean pi n and variance -# 2 , where 
/i n = E [Jj] = (n - l)dg n , 

< = Var [/J < (n - l)(2a^ n - (d ?n ) 2 ) < (n - l)(2<mg n ). 
Using the Central Limit Theorem [37, p. 183], we obtain 

Mi - /J, n \ <1p n } « 1 - Q 



(a) 

> 1 



7^ 



(A-4) 
(A-5) 

(A-6) 
(A-7) 



for all ip n > such that -0 n = o (n^firij- In the above equation, the Q(.) function is defined as Q(x) = 
^= / x °° e~ u2 l 2 du, and (a) follows from the fact that Q(x) < e~^, Wx > 0. Selecting ip n = (nq n ) 1 \/2t? n , 
we obtain 

P^-^Wl > 1-e-^. (A-8) 
Therefore, defining e = ^ = O ^(ng n )~l^, we have 

F{^ n (l-s)<I i < f i n (l + e)} > l-e~^\ (A-9) 

Noting that ng n — > oo, it follows that Jj ~ /x n , with probability one. Now, we show a stronger statement, 
which is, the contribution of the realizations in which | J 4 — (j, n \ > i\) n in the network average sum-rate is 
negligible. For this purpose, we give a lower-bound and an upper-bound for the network average sum-rate 
and show that these bounds converge to each other in the strong interference regime, when nq n — > oo. A 
lower-bound denoted by R^l, can be given by 



nWE 



> nWE 




Pihi 



h + 



n w 

M 

Pihi 



F{\Ii-fj, n \ <ij n } 



Ul + e) + ^f 



(A-10) 
(A-ll) 
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which scales as ^ log n (as shown in the proof of Theorem 1, by optimizing the power allocation function). 
An upper-bound for the network average sum-rate, denoted by R^l, can be given as 



nWE 



log ^1 
log (l 

< jfi + nWE 

(6) 



M 



T , NpW 

* ' M 

lo § I 1 + iw 



3 -(ng n )i 



k - Vn\ > ipn} 



N W 
M 



-(nq n )^ 



R^l + WO{nq n \ogn)e- {n<1 ^ 
& R^ 



(A-12) 
(A-13) 

(A-14) 

(A-15) 
(A-16) 



In the above equation, (a) follows from the fact that log(l + x) < x, (b) comes from the facts that 
^{Pihu} ^ Qn^ogn (this is shown in the proof of Theorem H]) and is fixed, and finally, (c) results 
from the fact that as nq n — > oo, nq n e~^ nqn ^ — ► 0. The above equation implies that substituting lj by its 
mean {{n — l)aq n ) does not affect the analysis of the network average sum-rate in the asymptotic case 
of K — > oo. 

Appendix II 
Proof of Lemma [2] 

Denoting Tj = {I E Cj \ hu > hrh}, the cardinality of the set Tj is a binomial random variable with 
the mean nF{h a > h Th }. From ©, we have 

M 



Rave ^ ^ 
3=1 



(B-l) 



where 



E 



5>(p w ,/f) 



zee,- 



E 



+ E 



zeT? 



(B-2) 
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in which denotes the complement of T,-. Note that 



E 



leJi 



n—E 

M 



M 



log 1 + 



log 1 + 



hupi 



T I N W 

11 T M 



hi > h Th 



W{hu > h Th } (B-3) 



hi, 



n w 

M 



hi > h 



Th 



V{h u > h Th ] 



< —e~ hTh E [h a \h u > h Th ] 
iv 



n 



—e-^(l + h Th ) 
i\ 



(B-4) 

(B-5) 
(B-6) 



where (a) follows from log(l + x) < x, for x > 0. It is observed that for hxh = clog n, where c > 1, the 
right hand side of (|B-6I) tends to zero as n — > oo. Thus, 



lim E 

n— >oo 



leJi 



0. 



Consequently, 



M 

lim VE 



and this completes the proof of the lemma. 



(B-7) 



(B-8) 



Appendix III 
Proof of Lemma [3] 



Using (1321) . we have 

E[Ei(pi, ha)] 



W 
MA 



^[huPi] + n 



MX' 



6rj 
X* 



nil — a 



'ma 7 



(a) W 



MX 
naW6r] 

MA' 4 
(ft) 14/ 



1 + r n )g„ 



naW . x 2 naW2n 



MX' 2 



MA' 3 



1 + Otfn 



Md 



\ 2 /-, \ 

.1 + T n )q n + Jfyi 1 + T n)Qn 

6/ 



6 

ra 2 



1 + r n + ^(1 + r n )e r " - ^(1 + r n )e 2r " 



2k 



(C-l) 

(C-2) 
(C-3) 



where £i = — - and £2 — — ttt- In the above equation, (a) follows from the fact that ElhupA = // 
roa war 



(1 + r n )q n , and (6) results from i) A = (n — l)aq n + 



n w 

M 



naq n and A' ~ naq n incurred by the fact that 



A > 1, and ii) q n = e~ Tn . Since naq n — > 00, it follows that the right hand side of (|C-3I) is a monotonically 
increasing function of r n , which attains its maximum when r n takes its maximum feasible value. The 
maximum feasible value of r n , denoted as f n , can be obtained as 



nae 



00 



logn. 



(C-4) 
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w 

Thus, the maximum achievable value for E[Hi(pi, hu)] scales as -j^-z logn. 

Appendix IV 
Proof of Lemma [4] 

i) Using © and assuming that all users follow the on-off power allocation policy, E[^(pj, ha)] can be 
expressed as 



leCi 



E[u l (p t ,h tt )] = J2^\Ri(P U \c ilh 



j = l,...,M, 



(D-l) 



where the expectation is computed with respect to hu and 1\. Noting that q n = ¥ {hu > r n }, we have 
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^(P (i) ,A (i) ) 



E 
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q n E 



^(P a) ,A (i) : 
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/iH < r n 
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+ (l-g n )E U(P (j) ,A (i) ) 



(D-2) 



hu < T n 



Since for hu < r n , pi =0, it is concluded that 



E 



M 



E 



log 1 + 



M 



For large values of K, we can apply Lemma \T\ to obtain 



E 



M 

QnW 
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where the expectation is computed with respect to hu- Using the Taylor series for log(l + x), (|D-5I) can 
be written as 
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J2-rbr- E [K\hu>r n ] 



M ^ k\ k 
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(a) q n W 
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^ k(naq n 
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where (a) follows from the fact that for large values of n, A « naq n . Also, (b) results from the fact that 
under a Rayleigh fading channel model, 

E [h u \h u > r n ] = 1+T n , (D-ll) 



2S 



E [h k u | h a > rj = r k n + kE [/iff 1 \ h u > r n ] . (D-12) 

Since A > 1, the term E [ fe » j^ a>Tn ^ <C E [ h u \^ l>Tn ] ^ w hi cn implies that we can neglect this term and 

simply write E [/i^l hu > r n ] ~ r*. (c) results from g n = e _Tn . Thus, (ID- II) can be simplified as 

ne~ Tn W ( t e Tn \ 

E[w<(pi, ^)] « u log 1 + • (D-13) 
M \ not J 

In order to find the optimum threshold value: 

f n = arg max E^^pij/iij)], (D-14) 
we set the derivative of the right hand side of (|D-13I) with respect to r n to zero: 

- e-- log (l + ^C-) + } + fr \ = 0, (D-15) 

\ not J na + r n e Tn 

which after some manipulations yields 

f n = logn - 21oglogn + 0(1). (D-16) 
ii) Using (ID- 161) . it is concluded that 

q n = e~ Tn (D-17) 
log 2 n 

= 5^-11, (D-18) 



n 



where 5 is a constant. 

Hi) Using (ID-161 ), we have 



e - — , (D-i9) 



na \ log n 

which implies that the right hand side of (ID- 131) can be written as 

EH (EES) « 2^. (D-20) 

Thus, the maximum value for E[ui(pi, hu)} in (|D-13I) scales as — — ■ logn. 

Ma 

Appendix V 
Proof of Theorem [3] 

Let us define Aj as the set of active links in cluster j. The random variable nij denotes the cardinality of 
the set Aj. Noting that for M = Q(K), lirn^^oo ^ is constant, it is concluded that n and rrij G [1, n] do 
not grow with K. To obtain the network average sum-rate, we assume that among M clusters, Y clusters 
have nij = 1 and the rest have nij > 1. We first obtain an upper bound on the average sum-rate in each 
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cluster when rrij — 1, 1 < j < M. Clearly, since only one user in each cluster activates its transmitter, 
Ii = 0. Thus, by using (|48l) , the maximum achievable average sum-rate of cluster Cj is computed as 



m = — e 

ave M 



, , ■ M r 



(E-l) 



where /i max = max {hu} ieC is a random variable. Since logx is a concave function of x, an upper bound 
of (|E-1I) is obtained through Jensen's inequality, E [logx] < log(E [x]), x > 0. Thus, 



(E-2) 



Under a Rayleigh fading channel model and noting that {ha} is a set of i.i.d. random variables over 



i G C,-, we have 



^(y) = P{/i max <2/}, ?/>0 

= l[nhii<y} 

ieCj 

= (i-e~y) n , 



where F h 



is the cumulative distribution function (CDF) of /i max - Hence, 

,n-l 



E[/i E 



o 



m/e y (l — e v ) dy. 



Since (1 — e 1 < 1, we arrive at the following inequality 



E[h v 



< 



nye v dy = n. 



Consequently, the upper bound of (|E-2I) can be simplified as 
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For rrij > 1 and due to the shadowing effect with parameters (a,w), the average sum-rate of cluster 



Co can be written as 



m = t -e 



log 1 



M 



(E-9) 



where i^'s are Bernoulli random variables with parameter a. Thus, 
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where is the sum of I i.i.d random variables {Z,i}\ =1 , where Zj = f3 ki h ki , k ^ i. For nij > 1, is 
greater than the interference term caused by one interfering link. Thus, an upper bound on the average 
sum-rate of cluster C, is computed as 



log 1 



F 

~NqW 

M 



+ 



3 — - 

YL v V i mj ~~ 1 ^ ; 
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log 1 + 
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where F = /i max = max {/iji} l6C .. According to binomial formula, we have 
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We have 
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Defining Z = (3 min hki and X = ^, the CDF of X can be evaluated as 

= ¥{X<x}, x>0 
= P{F < Zx} 

P{F < = z}f z (z)dz 



= J 

Jo 



fOO 



1 -- , 

e ^min 0(2; 
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= / (l-e-^fe-Ut. 
Jo 

Thus, the probability density function of X can be written as 

dFx (x) 
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Using the above equation, the right hand side of (IE- 151) can be upper-bounded as 

Y 
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log 1 



< n(3 r 



fx(x) log(l + x)dx 
"°° log(l + a; A 
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where the last line follows from the fact that < [3 min < 1. Substituting the above equation in (|E-14|) 
yields 



r& < ^(i-«r rl E 
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where (a) follows from (IE- 8 1) and the fact that m, e {2, n} does not scale with AT. 

Let us assume that among M clusters, T clusters have rrij = 1 and for the M — T of the rest, the 
number of active links in each cluster is greater than one. By using (IE-81) and (IE-291) , an upper bound on 
the network average sum-rate is obtained as 



R n 
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log 1 + 
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To compare this upper-bounded with the computed network average sum-rate in the case of M = 1, we 
note that as w < 1 and a < 1, we have a < 1 and consequently, 
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log 1 
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To prove that the maximum network average sum-rate obtained in (IE-321) is less than that value obtained 
for M = 1 from (fT8~l). it is sufficient to show 
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Since a < a, it is sufficient to show that rrij(l — a) mj 1 < -. Defining A(a) = am 3 (l — a) mj 1 , we 
have 

^ = m j (l-a) ffl '-- 2 (l-am j ). (E-36) 
da 

Thus, the extremum points of A(a) are located at a = 1 and a = where rrij G {2, n}. It is observed 
that 

A(l) = < 1, (E-37) 

and 

A (i)j^ir<, (E . 38) 

Since A(a) < 1, we conclude (IE-341) . which implies that the maximum average sum-rate of the network 
for M = Q(K) is less than that of M = 1. Knowing the fact that for M = o(K), similar to the result 
of Theorem 1, one can show that the maximum average sum-rate of the network is achieved at M — 1, 
it is concluded that using the on-off allocation scheme, the maximum average sum-rate of the network is 
achieved at M = 1, for all values of 1 < M < K. 

Appendix VI 
Proof of Remark 4 

Using © and © and for every value of 1 < M < K and a = 0, the average sum-rate of the network 
is simplified as 

\w ( hu : W 

(F-l) 
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Rave = ^2 ^ 
j=l ieCj 



w K 

\ M 



where the expectation is computed with respect to ha. Under a Rayleigh fading channel condition and 

K_ 

M ' 



using the fact that n = j-, (IF- II) can be written as 



f°° ( M \ 
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-e m Ei ~~~ (F-3) 
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e — / & (F-4) 

M Jit 

where Ei(x) = — Ei(— x) = s -j-dt, x > 0. Taking the first-order derivative of (IF-41) in terms of M 
yields, 

dR n „ P KW n$w_ ( N Q W\ [N W\ KW 



dM M 2 V M J \ M J M 2 
dR 

Since for every value of N W, is negative, it is concluded that the network average sum-rate is a 

monotonically decreasing function of M. Consequently, the maximum average sum-rate of the network 
for a = and every value of 1 < M < K is achieved at M — 1. 
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Appendix VII 
Proof of Remark 5 

From © and ©, the average sum-rate of the network is given by 
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where the expectation is computed with respect to ha. Under a Rayleigh fading channel condition, we 
have 
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To simplify (|G-4I) . we use the following series representation for Ei(x), 

_iy+i x s 



E 1 (x) = -7 + log C~)+jt 

\ / c 1 



x > 0, 



S.S! 



where 7 is Euler's constant and is defined by the limit [34] 
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0.577215665... 



Thus, (IG-41) can be simplified as 
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Consequently, the network average sum-rate for M = K is asymptotically obtained by 
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